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Abstract 

Using dynamic and thermodynamic magnetization measurements, we analyze the phase diagram of 
Bi2Sr2CaCu2 08+a single crystals containing amorphous columnar defects created by heavy-ion irradia- 
tion. Reversible magnetization experiments yield the respective magnitudes of the pinning energy and 
entropy contributions to the free energy of the vortex lattice. It appears that the entropy contribution in 
the London regime is relatively minor in both unirradiated and irradiated crystals, except in the case of 
high density of columns and inductions B that are smaller than the interaction field Hint ~ B^/6. The 
dependence of the entropy contribution on vortex- and defect density correlates well with measurements 
of the irreversibility line Hirr{T), which shows a sharp increase at Hint. 



1. Introduction 

The enormous enhancement of the critical current den- 
sity [1-4] and of the irreversibility line (IRL) separating 
the vortex glass and liquid phases [3,4] obtained after 
the introduction of amorphous columnar defects into 
the layered high-Tc superconductor Bi2Sr2CaCu208-i-i 
has spurred numerous investigations. Whereas initial 
studies conducted at temperatures below 20 K seemed 
to indicate that flux pinning and dynamics in this sytem 
could be described within a simple model describing the 
interaction between single pancake vortices (i.e. the in- 
tersections of the vortex lines with the superconduct- 
ing Cu02-double layers) and the insulating defects [1,2], 
later work showed that collective effects should be taken 



into account in order to understand the vortex physics 
at higher temperatures (above 40 K). The field depen- 
dence of the sustainable current density [5] can only be 
explained if the in-plane repulsion between pancakes is 
taken into account. Further experiments on the depen- 
dence of the screening current j on the angle between 
the applied magnetic field and the linear tracks demon- 
strated that above 40 K pancakes belonging to the same 
vortex line ("stack") move in a correlated way [6,7]; the 
stretched-exponential time dependence of the screening 
current, j ~ (Ini)^^''' [4,8], supports this finding. Mea- 
surements of the c-axis resistivity [9-11] and c-axis crit- 
ical current, using Josephson Plasma Resonance (JPR) 
[12-14], show a drastic enhancement of these quantities 
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when the field is increased above 5^/6; this has been 
interpreted in terms of pancake vortex aUgnment and 
an ensuing "recoupUng transition" in the vortex hquid 
phase. It was proposed that the recouphng is driven by 
the competition between the attractive interaction be- 
tween pancakes in difi^erent layers and the entropy gain 
obtained by spreading pancakes belonging to the same 
stack over difi'erent columns [14-16]. It must be stressed 
however that the fields under study are (much) larger 
than the crossover field Bcr = ^o/'y^s'^ (with 7 ~ 300 
the anisotropy parameter, s — 15 nm the spacing be- 
tween Cu02 double layers in the Bi2Sr2CaCu2 08+5 ma- 
terial, and $0 = h/2e the fiux quantum). This means 
that the interlayer pancake interaction is much weaker 
than the intralayer interaction. A number of authors 
have stressed the importance of the intraplane pancake 
repulsion in determining the vortex configuration and 
vortex dynamics in the presence of columnar defects 
[17,18]. In this paper, we analyze the reversible magne- 
tization of heavy-ion irradiated Bi2Sr2CaCu208-i-i sin- 
gle crystals. We determine the relative magnitudes of 
the free energy gain associated with vortex pinning on 
the irradiation-induced amorphous columnar tracks and 
the entropy gain associated with the possibility of pan- 
cake vortices to occupy different tracks. We show how 
information on vortex alignment can be obtained from 
the reversible magnetization, and how the vortex ar- 
rangement affects their dynamics and the phase dia- 
gram. 

2. Experimental details 

Different batches of Bi2Sr2CaCu2 08+a single crystals 
were grown using the travelling-solvent floating-zone 
technique at the University of Amsterdam and the Uni- 
versity of Tokyo. A number of crystals were annealed in 
air at 800C in order to obtain the optimum Tc = 90 K. 
Other crystals were left as-grown (slightly overdoped) 
with a Tc of 83 K. After imaging of the flux penetration 
by magneto-optics and the removal of defective regions, 
a number of large pieces of size 2x1 mm^ x20 mm 
were cut for the reversible magnetization experiments. 
Smaller pieces of size 800 x 800 x 20 nm^ were prepared 
for AC transmittivity measurements of the irreversibil- 
ity line (IRL) [8] . The crystals were subsequently irradi- 
ated with 5.8 GeV Pb^''^ ions at the Grand Accelerateur 
National d'lons Lourds (GANIL) in Caen. The ion 
beam was directed along the c-axis; each ion impact cre- 
ated an amorphous columnar track of radius cq ~ 3.5 
nm, traversing the entire crystal thickness. The ion dose 
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Figure 1: Magnetization of an optimally doped 
Bi2Sr2CaCu208-i-,5 single crystal grown at the Univer- 
sity of Amsterdam, and irradiated with 5.8 GeV Pb 
ions so as to produce a matching field B^ = 2 T. The ar- 
rows denote the irreversibility field Hirr as determined 
from the high temperature onset of a third harmonic 
response in field-cooled AC transmittivity experiments, 
as well as the interaction field Hint above which inter- 
vortex repulsion starts to limit the free energy gain that 
can be obtained from vortex localization on the colum- 
nar defects. 



for different crystals varied between 1 x 10^ and 2 x 10^^ 
cm^^, corresponding to dose-equivalent matching fields 
B^ = ^oUd between 20 mT and 2 T (n^ is the colum- 
nar defect density). The reversible magnetization was 
measured in Leiden using a commercial SQUID magne- 
tometer. The IRL's presented below 

were determined by ac Local Hall Probe Magnetom- 
etry [8,19] using an ac field of amplitude 6 Oe and fre- 
quency 7.753 Hz, applied colinearly with dc field and 
the crystalline c-axis. The irreversibility temperature 
Tirr{B) (or irreversibility fleld Birr{T)) is defined as 
that at which one first observes a third harmonic in the 
ac screening current when cooling. 

3. Results and Discussion 

S.l Reversible magnetization 

Figure 1 shows the magnetization of an optimally doped 
Bi2Sr2CaCu2 08-1-5 single crystal with B^ = 2 T. The 
reversible magnetization Mrev exhibits the usual fea- 
tures found after heavy-ion irradiation [20,21]. Note 
that the reversible magnetization corresponds to the 
vortex chemical potential /i, i.e. Mrev ~ —dG/dB = 
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— $0 — IdG/dUv = —^Q^n and is therefore equal to the 
free energy AG needed to insert (or to extract) a vor- 
tex into (from) the vortex lattice. At inductions much 
less than B,/, all vortices can lower their free energy by 
becoming localized on a defect site, hence the absolute 
value of the magnetization IMreul is reduced with re- 
spect to that of the usual magnetization in the London 
model, M°gj, = — eo/2'I'o \n{riBc2/eB), by a correspond- 
ing amount (eo(r) = $o2/47r^oA (T) is the vortex line 
energy, A(T) is the penetration depth, and rj « 1). The 
fact that this amount is large, i.e. of the same order of 
magnitude as the total magnetization, means that vor- 
tices are strongly bound to the columnar defects. As 
the matching field is approached, favourable sites be- 
come rare and more and more "interstitial" vortices {i.e. 
not localized on a columnar defect) with higher free en- 
ergy appear. Hence |Mrej;| increases (at the interaction 
field, labelled Hi„t in Fig.l). At fields far above B^, 
no defect sites are available anymore, all vortices that 
enter the sample occupy sites in the intercolumn space; 
hence, the magnetization is comparable to M^^^. The 
rounding of the magnetization near B^ is the result of 
intervortex repulsion, which prohibits the occupation of 
close-lying columnar defect sites even at fields well be- 
low B^ [18,20], and forces vortices to occupy interstitial 
sites. 

The free energy gain manifest at low fields was ini- 
tially interpreted as being the sole result of pancake vor- 
tex pinning on the columnar tracks. Then, Mrev{B « 
B^) — M°(.„ -I- f/o/$o, and the pinning energy per unit 
length C/o can be directly obtained as the difference be- 
tween Mrev (at low field) and the extrapolation to low 
fields of the magnetization at fields much above B^. 
The application of this procedure to irradiated crys- 
tals with relatively low B^ has been used to establish 
that the main pinning mechanism in heavy-ion irradi- 
ated Bi2Sr2CaCu2 0s+5 arises from the interaction of 
the vortex core with the defects [22] . Applying the same 
method to the magnetization of crystals with higher B^ 
(as in Fig.l) yields pinning energies that exceed any rea- 
sonable theoretical prediction. Conversely, in numerical 
simulations of column occupation [18] one has to assume 
an artificially short-ranged vortex interaction in order 
to arrive at a pronounced magnetization feature such 
as that in Fig.l. Furthermore, the fact that the loga- 
rithmic field derivative of the magnetization in the low 
field limit is substantially smaller than that at fields 
larger than B^ indicates that there must be a second 
free energy contribution (i.e. entropy) that favours vor- 
tex creation and column occupation. 
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Figure 2: Temperature dependence of the logarithmic 
field derivative of the reversible magnetization for the 
same crystal as in Figure 1, in the field regimes B <^ B^ 
and B 3> B^ respectively. The drawn lines represent the 
value of dM^^^/dlnB — eo(0)/2<J>o as extracted from 
experiment, and the hypothetical entropy contribution 
dMth/dTdlnB = fcsT/$o [23,24]. 



S.2 Entropy contribution to the magnetization 

The possible entropy contribution Mth to the magneti- 
zation of layered superconductors in the London regime, 
arising from vortex positional fluctuations, was consid- 
ered by Bulaevskii et al. [23]. They derived that for 
large fields B > Bar, Mth = (fcsT/$os) ln(Bo/B), 
with Bo ~ Bc2{£os/kBT) a characteristic field related 
to the elemental phase area. As a consequence, the 
logarithmic field derivative dM^^^/dluB is cancelled 
by dMth/dlnB at the "crossing temperature" T* — 
Tc/{1 + 2kBT/eos); at T* , the reversible magnetization 
is field-independent. The theory was extended [24] to 
describe heavy-ion irradiated layered superconductors, 
which exhibit two distinct crossing temperatures Ti and 
Ta* [16,21], in the field regimes B « B^ and B » B^ 
respectively (see Fig. 2). 

The magnitude of the entropy contribution 
to the magnetization can be unambiguously 
evaluated from the temperature derivative of 
dMrev/d\nB. Namely, this is expected to con- 
sist of the two contributions d^Mrev/dTdlnB = 
-eo(0)/2$oTc and d^'Mth/dTdlnB = ke/^os 
[25]. The first term can be read from Fig. 2, 
d^M°^^/dTdlnB ^ 23, and yields the value of 
the penetration depth at T = 0, A(0) = 180 nm [26]. 
The second contribution is constant and equal to 4.4 
Am^^K^^; the ratio 4.4/23 = 0.19 represents the error 
in A(0) due to the presence of fiuctuations. Having 
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Figure 3: (a) Reversible magnetization of a series of Ughtly 
overdoped single crystals (University of Amsterdam, Tc = 83 
K), irradiated to varying matching fields, as indicated. The 
drawn lines indicate the magnetization of the unirradiated 
sample, and fits to the data for B^ = 0.2 T and 0.5 T to 
Eq. (1) (b) The difference between the experimental magne- 
tization and that of an unirradiated crystal, multiplied by 
#0*- In the low-field limit, this quantity is equal to the sum 
of the pinning energy and the entropy per pancake. The 
drawn lines indicate the "bare" pinning energy per pancake 
UOs (line through data for B^ = 0.5 T), and, for higher B^ , 
the sum of Uos and the entropy TS = (fcsT/<I>os) \n{/3B^/B) 
associated with the fact that individual pancakes belonging 
to the same vortex line can take advantage of several colum- 
nar defect sites. The interaction field Hint is identified as 
that at which, as a result of intcrvortex repulsion, pancake 
vortices belonging to the same stack start to line up on the 
same columnar defect. 
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Figure 4: Matching field dependence of the low-field en- 
tropy contribution to the free energy (measured at an applied 
field fioHa = 0.1 T). The drawn lines indicate the model 
dependence in the regime where pancakes belonging to the 
same stack are ordered on the same column over a length La 
( S ^ B^\Ti{f3B^/ B)) and that where vortex lines are dis- 
sociated into single pancakes occupying different columnar 
defects (S ~ ln(/3B^/B)). The points for B^ = 0.2 and 0.5 
T are in fact upper limits for TS, obtained from the error 
bar in the magnetization measurements. 



established tlie value of A(0) with confidence, we 
use it to evaluate the theoretical logarithmic slope 
dM^^JdlnB = £o(r)/2$o (drawn in Fig. 2). The 
close coincidence with the data points for B S> B^ indi- 
cates that in this field regime the entropy contribution 
is smaller than the experimental error bar. The near- 
equality of Mrev{B ':^ B^) with the magnetization of 
unirradiated crystals indicates that in either case the 
entropy in the London regime is very modest indeed, 
and that vortex positional fluctuations account for only 
a small correction to the total free energy. 

The situation is very different for the regime B <^ B^, 
in which all vortices are localized on a columnar defect, 
and where the experimental dMrev/dlnB lies clearly 
below eo(T)/2$o. Although the difference is smaller 
than the value kBT/^os expected from theory [24], the 
fact that it exists indicates that entropy is important 
in the low fleld regime. This is further borne out by 
Fig. 3, which shows the reversible magnetization of a 
series of lightly overdoped crystals with different num- 
bers of sites available to vortices, i.e. with different 
defect density Ud or matching field B^. It appears that 
Mrev strongly depends on the density of pinning sites 
at low fields {B <^ B^), but not at high fields. A strong 
B^ -dependence of the low-field magnetization cannot 
be explained in terms of pinning only, but is consistent 
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with a substantial entropy associated with the possibil- 
ity for pancake vortices to occupy different columnar 
defect sites. 

The matching-field dependence of the entropy may 
be extracted as follows. First, subtract the magnetiza- 
tion of the unirradiated crystal from the experimental 
curves. In the limit B <^ B^ this yields a quantity, 
which, when multiplied by <^os, is exactly equal to the 
sum of the pinning energy and the entropy per pan- 
cake vortex. Next, we observe that for relatively small 
matching fields, B^ < 0.5 T, dMrev/dlnB in the low- 
field and high-field limits are nearly equal. This implies 
that for these matching fields, the entropy is small with 
respect to the pinning energy over the entire investi- 
gated field range. We can then confidently determine 
the pinning energy per pancake Uos, either by a fit to a 
suitable expression, e.g. that given in Ref. [20], 
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or by simple identification of the low- field plateau in the 
B^ = 0.5 T data of Fig. 3(b) with Uos. The entropy 
contribution to the free energy follows as the difference 
between the data and the pinning energy, i.e. TS = 
^os{Mrev - M°g„) - Uos. It is depicted in Fig. 4. We 
find that the entropy is small at modest B^, then rapidly 
rises between 5^ = 0.5 T and 1 T, after which there is 
a further increase proportional to InB^. Fig. 3 shows 
that in the latter regime TS is proportional to InB. 

We interpret this behaviour as follows. At very small 
matching fields B^ <^ Bcr, the (attractive) electro- 
magnetic and Josephson interaction between pancakes 
in adjacent layers prohibits vortex wandering between 
columnar defect sites; in this regime entropy can only 
be gained at the expense of pinning energy by vor- 
tex wandering into the intercolumn space [27,28]. Ev- 
idence for this entropic reduction of the pinning en- 
ergy for very small matching fields was obtained in 
Refs. [7,8]. As the density of sites is increased (to 
above 6Bcr, see below), vortex lines can gain free en- 
ergy by spreading to neighbouring columns. The en- 
tropy gain in this regime occurs at the expense of elas- 
tic (shear and tilt) energy only, and can be estimated 
as TS = ksTlnW = ksTlnNcLs/La, where Ls is the 
thickness of the sample in the field direction. La is the 
distance over which pancakes belonging to the same vor- 
tex line remain ordered on the same columnar track, and 
A'^c = I3B^/B is the number of columns available to any 
one vortex. One then has 



Mth 



kgT s ^^(PB^\ 



(2) 



From the extrapolation of the low field data in Fig. 3 (b) 
to Uos we obtain /3 ~ 0.5. Balancing the gain in thermal 
energy and the loss in tilt deformation energy gives La ~ 
eo/ifUdkBT) and 
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in agreement with the rapid rise between B^ = 0.5 T 
and B^ = 1 T (see Fig. 4). The typical stacklength 
can be obtained as the ratio between the experimental 
entropy values and {ksT /^qs). As the defect density 
increases. La ^ s{eos /kBT)[Bcr / B^) decreases, until 
at B^ ^ Bcr{eQs/kBT) it becomes equal to the layer 
spacing s. Experimentally this crossover occurs at a 
temperature dependent matching field. For T = 74 K, 
it is at B^ ~ 1 T. When the alignment length has de- 
creased to a single pancake vortex, further entropy gain 
can only be obtained through the supplementary defect 
sites made available to the pancakes at still higher B^; 
then, following Ref. [24], 



Mi, 






(4) 



in very good agreement with the experimental data 
shown in Figs. 3 and 4. 

The field dependence of the reversible magnetization 
is explained as follows. For the lower matching fields 
B^ < 0.5 T the entropy contribution is negligible with 
respect to the pinning energy, hence the increase of 
Mretil at Hint (labelled in Fig. 3) must be ascribed 
to the reduction of the total pinning energy due to in- 
traplane pancake repulsion. This can be accompanied 
by the appearance of the first interstitial vortices. At 
higher matching fields, the low-field magnetization is 
the sum of Mre„, the gain in pinning energy per pancake 
Uos, and the entropy gain (4) arising from the fact that 
each pancake can occupy different columns. As the field 
is increased, less sites are available per pancake, and the 
entropy contribution decreases logarithmically, until, at 
Hint, intraplane repulsion forces pancakes belonging to 
the same stack to lign up on the same column. This 
causes the rapid drop of Mth, which is now described 
by Eq. (2), and the increase of |Mrei,|. As the match- 
ing field is approached, the increase of \Mrev\ is further 
enhanced by the decrease in pinning energy and the ap- 
pearance of interstitial vortices. 
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Figure 5: Irreversibility line for an optimally doped 
Bi2Sr2CaCu2 0g single crystal (University of Amster- 
dam) irradiated with matching field B^ — 0.5 T. The 
IRL is determined as the temperature at which a third 
harmonic in the AC shielding current can be first ob- 
served upon cooling. The drawn line represents a fit to 
Eq. (5) with the indicated parameter values. 
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Figure 6: Irreversibility lines for optimally doped 
Bi2Sr2CaCu208 single crystals (University of 
Tokyo) with different matching fields < B^f, < 
T. The line labelled BpoT indicates the first order 
transition field in the pristine crystal. Other lines 
are guides to the eye. 



3.4 Relation with the phase diagram 

In the following we dicuss the consequences of the pan- 
cake vortex arrangement for the mixed state phase di- 
agram and vortex dynamics. Fig. 5 shows the IRL 
determined on an optimally doped Bi2Sr2CaCu208+i 
crystal with B^ = 0.5 T, together with the //int-line 
above which intraplane vortex repulsion limits the free 
energy gain due to columnar defect pinning. The Hint- 
Vine depicted in Fig. 5 agrees very well with the "recou- 
pling transition" line found in the JPR experiments of 
Refs. [f3,f4]. 

The irreversibility line has three well-defined parts. 
At low fields, HirriT) increases exponentially 
as temperature is lowered. When the interaction (or 
"recoupling" ) field Hint is reached, the exponential in- 
crease stops and Hirr rises sharply. Thus, at the precise 
field where intraplane vortex repulsion starts to deter- 
mine the vortex arrangement in the columnar defect 
potential, vortex localization is spectacularly improved 
and vortex dynamics slowed down. Above Hint, every 
vortex line is localized not only by the defect poten- 
tial but also by the potential created by the neighbour- 
ing vortices. As field is increased further, Hirr even- 
tually starts to decrease roughly linearly with temper- 
ature, presumably because the proliferation of weakly- 
localized interstitial vortices determines flux dynamics. 

The behaviour of the IRL for different matching 



fields is shown in Fig. 6. In the low-field, "exponen- 
tial" regime, Hirr increases monotonically with ion dose 
up to B^ ~ 30 mT, after which it saturates to a 
dose-independent exponential dependence on temper- 
ature. Such a saturation of the irreversibility field 
is reminiscent of experiments on heavy-ion irradiated 
YBa2Cu307_f , where it signals the stability limit of the 
entangled vortex liquid with respect to the introduction 
of correlated disorder. In the low-matching field regime 
below the saturation (i.e. B^ < 30 mT), the exponen- 
tial temperature dependence of Hirr can be observed 
up to Hirr B^, after which it crosses over to a linear 
dependence Hirr ~ (1 - T/Tc). For B^, > 30 mT, the 
abovementioned break at Hint ~ B^/6 appears. It is 
found for the whole investigated range 0.1 T < B,/, < 4 
T, always near the same fraction B^/6. At fields above 
B^/6, Hirr decreases with increasing matching field. 
This is 

again the result of the fact that an increase of the num- 
ber of available columnar defect sites promotes vortex 
wandering and eventually, diffusion. 

3.5 Vortex delocalization in the low- field limit 

The low-field regime H < Hint in which the irreversibil- 
ity line saturates to a B^- independent exponential tem- 
perature dependence is that in which interactions be- 
tween different vortex lines are irrelevant, and, there- 
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fore, where every pancake vortex can find a suitable 
columnar defect site. Moreover, since for B^ < 0.5 T, 
pancakes belonging to the same vortex line are well- 
aligned even in the vortex liquid {La ~ 10s at 74 
K), the coincidence of the exponential part of the ir- 
reversibility lines for 0.1 T< B^ < 4 T implies that 
near HirriT) pancakes belonging to the same stack are 
aligned on the same ion track irrespective of the match- 
ing field. Vortex motion and delocalization can then 
be modelled in terms of a "discrete superconductor" 
in which the only sites allowed to the vortices are de- 
fect sites. The absence of strong positional fluctuations 
inferred in Sections 3.1 and 3.2 implies that the main 
thermal fluctuations will be occasional single-pancake 
jumps to neighbouring columns. These can be seen 
as bound vacancy- interstitial pairs (or "quartets" [30]), 
with separation R and energy £q = £os{R/'ys)^ , present 
in the vortex system of the discrete superconductor. 
The onset of vortex diffusion at the IRL corresponds 
to the defect pair unbinding; hence, ksTirr ~ eq{Ri) 
where Ri ~ (<l>o/-B)^^^ exp(-£os/2fcsr) is the typical 
distance between free dislocations (free pancakes) in 
the vortex liquid. Gathering terms, one has kaTirr = 
eos($o/7^s^)exp(eos/fcsr) or 



JDir 



Be 



SqS 

ksT 



exp 



(S)<«- 



<C -B < Hint 



(5) 



which satisfactorily describes the exponential part of 
the IRL, with parameters A(0) = 165 nm and Bcr ~ 
6.5mT (i.e. 7 ~ 350) (see Fig. 5). The upper limit of 
applicability of Eq. (5), as this emerges from Section 
3.3 {i.e. B < Hint = B^/Q) also gives the lowest B^ 
-value for which Eq. (5) applies: B^ = QBcr ~ 40 mT, 
in very good agreement with the matching field value 
at which the IRL saturates. We conclude that the value 
B^ = 6-Bcr delimits the low defect-density regime, in 
which thermal vortex excursions into the intercolumn 
space are allowed, from the high-defect density regime 
where the "discrete superconductor" -model of strongly 
localized vortices is applicable. 

4. Conclusions 

Reversible magnetization measurements on 
Bi2Sr2CaCu2 08+,s single crystals containing amor- 
phous columnar defects introduced by heavy-ion 
irradiation are used to determine the magnitude of 
the pinning energy and entropy contributions to the 
mixed state free energy. Whereas the pinning energy 



contribution can always be identified, the importance 
of entropy depends strongly on the defect density 
and the magnetic field. For low defect densities and 
unirradiated crystals, the entropy contribution to the 
free energy in the London regime {i.e. the entropy 
gain associated with vortex positional fluctuations) is 
found to be very minor with respect to the vortex core, 
electromagnetic, and pinning energy contributions. It 
therefore does not determine the main features of the 
phase diagram. At higher matching fields (typically 
above 1 T) and low magnetic fields B < Hint B^/6, 
the entropy associated with the possibility of pancakes 
belonging to the same vortex line to occupy different 
defect sites becomes important. This configurational 
entropy greatly exceeds the fiuctuation entropy in 
unirradiated samples; it can attain values that are 
more than twice the pinning energy, for B^ ~ 4 T. 
When the magnetic field is increased above Hint, 
intraplane pancake vortex repulsion causes the total 
pinning energy and the configurational entropy to 
diminish. Nevertheless, vortices become more strongly 
localized, since they are now retained not only by 
the columnar defect potential, but also by that of the 
neighbouring vortices. This effect is apparent as a 
steep increase in the IRL, and as the recently found 
"recoupling" transition in JPR experiments [13,14]. 
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